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Abstract 

In the Littlest Higgs model with T-parity (LHT), top-partners (T_) are produced in pairs at the 
Large Hadron Collider (LHC). Each top-partner decays into a top quark (t) and the lightest T-odd 
gauge partner Ah- We demonstrate reconstruction of the ti system decaying hadronically, and 
measurement of the top-partner mass from the tot2 distribution. A top quark from a T_ decay 
is polarized, and we discuss the effect of this polarization on the decay distributions. Because the 
events consist of highly collinear jets which occasionally overlap each other, we compare distribu- 
tions using different jet reconstruction algorithms (Snowmass cone, kt, Cambridge, SISCone). We 
find clustering algorithms are advantageous for studying top polarization effects. 
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I. INTRODUCTION 



In spite of the success of the Standard Model (SM) in explaining the particle interactions, 
there are yet two unsolved questions in the SM. One is the fine tuning problem and the other 
is the existence of the dark matter (DM) in our universe. 

The fine tuning problem is the question why the Higgs is likely to be so light as expected 
from the LEP data (m^ < 198 GeV at the 95% confidence level while naturally the 
Higgs mass is the order of the cut off scale of the theory due to radiative corrections. Some 
mechanism should protect the Higgs mass, and interactions involving the Higgs sector should 
be extended from the SM one. 

The existence of the DM is now^ established by the cosmological observations such as 

2. 



WMAP, SDSS and SN-Ia 



. The SM should be extended to include the DM, 
which should be a neutral stable particle. Moreover, provided that the DM is a thermal relic 
whose strength of the coupling to the SM particles of the order of the weak interaction, the 
mass can be a few 100 GeV from a rough dimensional analysis. 

To solve the fine tuning problem, many models have been proposed. However, a so-called 



'little hierarchy problem" 



aa 



arises in these models, once constraints from precision 



measurements are imposed. New operators arising from these models naively have the cut 



off scale A > 10 TeV to be consistent to the experimental data 



9]. 



Some successful models solving the little hierarchy problem have parity structures, for 



example t 



m, 



14 



re MSSM with R-parity 



15 



16 



ity Q 
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l^ . the little Higgs model with T-parity (LHT) 



171 ]. or the universal extra dimensional model (UED) with KK-parity la. Il9|. 



Such a model which has a parity structure predicts the stable lightest parity odd particle as 
a candidate of the DM. The cut off scale of new operators in the model can be as low as a 
few TeV [2o| , and the mass of the lightest parity odd particle can be on the order of a few 
100 GeV. 

The LHC is starting soon, and it is likely to discover new particles with masses up to a 



few TeV 



21 



22]. The models mentioned above are studied intensively by many authors. 



All of these models have partners of the SM colored particles which decay into the stable 
lightest parity odd particle through the parity conserving interactions. 

At the LHC, these partners are produced in pairs. The signal is multiple high p? jets and 
high pt leptons, each accompanied by missing transverse momentum J^t- Among the signals 
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studied so far, the signals with high pt leptons are very promising 23j, |24j , because the SM 
backgrounds are smaller. However, branching ratios of new particles into leptons strongly 
depend on model parameters. In addition, Jets + leptons signals are often accompanied by 
undetectable neutrinos, which are also the source of missing momentum. They sometimes 
reduce the significance of the kinematical endpoints of the signal distributions for mass 
reconstructions. And even in the case that the lepton branching ratios are small, events 
with multiple jets and no lepton are enormously produced. 

In this paper, we focus on the top partner signal in the LHT. The top partner pair 

n 

productions occur with sizable cross section [25j compared with stop pair productions in 
the MSSM because the top partner is a fermion. The event has simple kinematics, ti + Ifi? 
where top quarks are highly boosted 1 . Decay products from a boosted top are collimated 
and they are easily identified as a jet system with mass ~ m t with high probability using 



hemisphere analysis 



27 



281 ]. On the other hand, by imposing lepton veto and top tagging, 



background events from tt+jets can be reduced significantly. 

The process above was partly studied in Ref. 28|. In the paper, COMPHEP 29] and 
HERWIG6 . 5 {30]] are used for the event generation and AcerDETl . 32I is used for the de- 
tector simulation and jet reconstruction. AcerDETl . implements the Snowmass cone al- 
gorithm. However, the algorithm is not optimized in resolving the overlapping jets arising 
from boosted top quark correctly. 



Recently, Fast Jet 



34 



35 



36 



331 ] was re 



37J, Cambridge 



38 



eased. Infrared stable jet reconstruction algorithms (kt 



391 ] . SISCone |40j]) are implemented in the code, which 



is significantly faster than previous codes. We study improvements with these advanced 
jet reconstruction algorithms. In this study, we interface AcerDET calorimeter information 
to Fast Jet and reanalyze the same process as in Ref. [28] to compare the results. We 
find the kt and Cambridge algorithm have advantage to resolve overlapping jets. We also 



generate the signal events with underlying events using HERWIG6.5 + Jimmy 41], and find 



the jet resolution with the kt algorithm is significantly affected by underlying events. This 
motivates us to show results with the Cambridge algorithm mainly in this paper. 

We study the potential to measure the mass of the top partner using the reconstructed 



1 Similar highly boosted tt signals have been discussed in the RSI model 26|. However, it is different from 
the signal considered in this paper, because it is not accompanied by I$t- 
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top candidates, although only the discovery potential is discussed in Ref. [281] . One of the 



important variable is a rriT2 



421 ] . The rriT2 is a function of two visible momenta, a missing 



transverse momentum and a test mass. In the case that the mass of the lightest T-odd 
particle (?tt.ltp) is known, the endpoint of the txiti distribution is equal to the top partner 
mass at M test = rrii^p. Therefore we can measure the top partner mass mrusing this 



43J+ HERWIG 



distribution. We also generate the Standard Model backgrounds using ALPGEN 
and conclude that they do not affect the endpoint of the rriT2 distribution. 

We also discuss top polarization effects. A typical LHT model predicts a top partner 
which decays dominantly into a right-handed top quark tn and a heavy photon Ah- The 
polarization of tops can be measured through investigating decay distributions of tops. And 
we show that there are distinguishable difference between completely polarized case and 
non-polarized case in jet level analysis. 

This paper is organized as follows. In Sec. [Til we explain our simulation setup for studying 
a top partner at the LHC And we show how to reconstruct momenta of top quarks arising 
from top partner decays and how to measure the top partner mass using a reconstructed 
m,T2 distribution. In Sec. IIIIl we discuss differences among the jet reconstruction algorithms. 
In Sec. HVt we study top polarization effects. Sec. [V] is devoted to the discussions and 
conclusions. 



II. TOP PARTNER RECONSTRUCTION AT THE LHC 



A. Event generation 

In the following, we assume the top partner is the lightest in the fermion partners and 
decays exclusively to the lightest T-odd particle Ah and a top. The top partner may be 
produced in pairs at the LHC and decays as, 

pp -> T_T_ -> tiAjjAjj -> bW + bW-A H A H -> 6j + tf T . (1) 

This process is similar to scalar top (t) pair production process in the MSSM. The production 
cross section of top partner is larger than that of scalar top in the case that the masses are 
the same, because top partner is a fermion. At the LHC, T_T_ production cross section is 
0.171 pb for rriT_ = 800 GeV. In order to identify this process, it is important to tag top 
quarks, and measure missing transverse momentum I$ T arising from escaping A#'s. 
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In the previous study 28J], the events pp — > ttAuAu are generated by COMPHEP 291, and 
top q oa rk _ are frL* to HERWIG6 . 5 Q. In this pape r , we a ls0 use COMPHEP 
and HERWIG6 . 5 for signal events generation. Our study is based on 8,550 signal events 
corresponding to J dtC = 50fb _1 . 

The Standard Model background to the signal comes from the production of QCD, ti+n 
jets, W + n jets, and Z + n jets events. In ATLAS study, it was shown that these four 
processes contribute to the background of the lepton +I$t+ jets channel for SUSY search 



with approximately the same order of magnitude [31]. Among those, the QCD background 
arises due to the detector smearing and inefficiency, and we do not attempt to simulate it 
in this paper. Even if QCD background is taken into account, it will not affect the results 
significantly because we require top mass cuts for the event selection. We will discuss this 
point later. The other processes contributes to the background due to hard v produced 
from Z and W decay. They are generated by ALPGEN+HERWIG in the paper. To reduce 
the computational time, we generate Z{— > vv) + n jets (n < 4) events corresponding to 
5 fb _1 with J$t > 150 GeV, W{— > Iv) + n jets (n < 4) corresponding to 10fb -1 with 
^partem &r > 400 GeV , and ti+n jets (n < 2) corresponding to 12 ftr 1 with ^partem > 
500 GeV respectively. Parton shower and matrix element matching are performed using 
MLM scheme provided by ALPGEN. We require ?7 max = 5, prmin > 30GeV and Rjj > 0.4 
for parton level event generation before the matching. The tiZ followed by Z — > vv events 
become irreducible backgrounds, and we have generated the events for roughly 50fb _1 . We 
do not apply K-factor both for signal and background. 

We use AcerDETl . for detector simulation, particle identification and jet reconstruc- 
tion as in [281 ] . In addition, we interfaced calorimeter information of AcerDETl. to 
Fast Jet2 . 2beta 33] so that we can compare different jet reconstruction algorithms si- 
multaneously, and model the detector granularity. Here the calorimeter information is the 
energy deposit ^ sum to each cell i centered at ((pi, r/i) with the size A(j) = 0.1 and Aw = 0.1 
(0.2 in the forward directions) without smearing. They are interfaced as massless particles 
with momenta pi = (E^ um ,r]i, (pi) to Fast Jet 2 . In this paper, we study jet distributions in 
the infrared stable algorithms, (kt, Cambridge and SISCone), together in those for the Snow- 
mass cone algorithm provided by AcerDET. To compare the four jet algorithms under the 



2 Effects of shower propagation to nearby cells are not taken into account. 
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same conditions in Sec. IIII A\ we switch off the jet energy smearing. For the background, the 
smearing on I$t might have the same impact on the estimation of the number of events after 
the cut, therefore we use the I$t smeared by the AcerDET. Effects of Jet energy smearing 
are discussed in the Appendix IA 1[ 

B. Event selection and Top reconstruction 
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TABLE I: Summary of the number of events after the cuts for JdtC = 50 fb . Numbers of BG 
events are properly scaled to 50 fb -1 . See the text for the detail. 



We describe our cuts to select T_T_ events. The summary of the numbers of the events 
after the cuts is shown in Table [H First, we impose our standard cuts for jet pt, $t and 
veto high p T isolated leptons, 

$ T > 200 GeV and tf T > 0.2M efr , n 50 > 4 and n 100 > 1, n lep = 0. (2) 

Here, 

M eff = E PT t0n + E PT° t0n + $T, (3) 

p T >50GcV p T >10GeV p T >10GeV 

t?| <3 l')l<2.5 \r/\<2.5 

n 50 (ft-ioo) is a number of jets whose pt is larger than 50 (100) GeV. n\ ep is a number of isolated 
leptons (e, fi) with pt > 5 GeV and \rj\ < 2.5. Missing transverse energy is calculated using 
the energy deposit to the calorimeter and isolated leptons. It is calculated with smearing 
for the Standard Model background calculation. 

The lepton cut reduces tt+ jets and W+ jets background, in which large 1$t is dominantly 
caused by neutrinos from leptonic W decay. W background still remains because W can 
decay into r. 
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We applied a hemisphere analysis to find top candidates [27J. Each of high jets 
(Pt > 30 GeV and < 3 ) is assigned to one of the two hemispheres which are defined as 
follows; 

Vz e Ht,j e H 2 d(p Hl ,Pi) < d(pH 2 ,Pi) and d(P H2 ,Pj) < d(P Hl ,Pj), (4) 

where 

P Hi = Pfc, (5) 

keHi 

\ (Ei-\Pi\ cos 9ij)Ei 
d{puVi) = {Ei - Ej)2 , (6) 

cos 6^- = . * . J " . (0^- is the angle between pj and pj). (7) 
iPillPil 

To find hemispheres, we first take the highest px jet momentum p± as Ph 1 and take the 
jet momentum pk which maximizes AR(pi,pk) ■ PkT among all k as pn 2 - We group jets 
into hemisphere Hi (i = 1,2) according to the eq.(Tj0). New P#'s are then calculated from 
Eq.((SD, and this process is repeated until the assignment converges. In this analysis, collinear 
objects tend to be assigned into the same hemisphere. Top quarks from T_ decays are highly 
boosted, then the decay products from the two top quarks are correctly grouped into different 
hemispheres with high probability. In this situation, the dependence on the definition of the 
distance eq.(jS]) is weak. Change the definition of the distance as d(pi,pj) = 9ij causes 
negligible differences of the acceptance and our analyses in this paper. 

To assure the correct top reconstructions, we require both hemispheres' transverse mo- 
menta are larger than a threshold, 

Pt,h x > Pt,h 2 > 20° GeV. (8) 
After imposing these cuts, distributions of the invariant masses of the hemisphere mo- 



menta (mp H = y Pff.) for the T_T_ and the Standard Model background events are shown 
in Fig. [U We can see peaks at the top mass both for the T_T_ and tt+jets events in mp H1 
distributions (Fig. [TK, anddb). On the other hand, such a peak is not seen in the mp H2 
distribution for tt+jets events (Fig. [TH). This is because at least one of the two tops should 
decay leptonicaly to give large J^t- We also plot the distribution for Z+ jets and W+ jets 



with dashed and dotted lines respectively. We do not see any structure in the hemisphere 
mass distributions. Two dimensional scattering plots in mp H] vs. mp H plane for the signal 
and tt+jets events are also shown in Fig. [2J which show the clear difference between them. 

We can reduce the events from Z+ jets, W+ jets and tt+ jets with hemisphere mass cuts. 
In Table HI the column mp H1 ~ m t shows the number of signal and background events after 
requiring 150 GeV < mp H1 < 190 GeV. The number of events of tt+ jets (Z+ jets, W+ 
jets) decreases by approximately 1/6 (1/13, 1/17) after the cut. In the second hemisphere, 
only the signal distribution (Fig. lb) has a peak and background distributions (Fig. Id) are 
fiat. 

We do not simulate QCD background in this paper. The magnitude of the QCD back- 
ground for SUSY 0-lepton channel is approximately the same as that of the tt+ jets back- 
ground jsi|. The hemisphere mass distribution of QCD background should be similar to 
that of Z+ jets or W+ jets background. Therefore the contribution of QCD background 
after the hemisphere mass cuts may be approximately the same as that of Z+ jets or W+ 
jets background. 

The column "both" shows the number of events after requiring the cut that both mp H1 
and mp H2 are consistent with m t . The number of background becomes small by a factor 
of ~ 1/100. However, the cut also reduces the signal events by a factor of 1/10. This is 
reasonable because we have the minimum jet energy cut for the hemisphere reconstruction 
(Pt > 30 GeV), and some of top decay products may not contribute to the hemisphere 
momentum. Additionally, in the case that a 6-quark decays semi-leptonically, the hemisphere 
momentum and the invariant mass are also reduced. Maybe the approach in Ref. 49] 
improves the mass resolution further. At this point, the background still dominate the signal, 
S/N ~ 1/3. If the sideband events can be used to estimate the background distribution, 
the significance of the signal events goes beyond 5 sigma. 

To verify whether a momentum of a hemisphere whose mass is consistent with m t correctly 
matches a top momentum, we compare the Ph 1 with momentum of top partons p pa rt- In 



Fig. El we show the distributions of the Apx = Pt,h x — Pr, P art and AR = a/ (Arj) 2 + (A0) 2 
to see the difference between the two momenta. Here, AR and Ap? is defined for one of the 
two top partons that gives smaller AR. The Apt/pt,H! distribution has peak near 0, and 
AR < 0.05 for most of the events. We conclude that a momentum of a hemisphere with 
rrip H ~ m t may be considered as a momentum of a top partons. 
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FIG. 1: The distributions of the invariant masses of a) H\ and b) H 2 for the T_T_ events for 50 
fb _1 , c) Hi and d) H 2 for the events from tt+jets (solid) , IV+jets (dotted) and Z+jets (dashed) 
after the cut Eq.© and Eq.©. BG events are rescaled to 50 fb _1 . 



C. Measurement of the end point of m,T2 

We now show that top partner mass (mr ) can be measured using the endpoint of the 
distribution of the Cambridge variable 42] for tt + I$t system if we know the LTP mass 
i m A H )- First, only the signal distribution is considered and after that we will show that the 
background event does not contribute to the events near the endpoint and can be neglected 
for the determination of the endpoint. It turns to be the best discrimination between the 
signal and SM backgrounds. 

This variable is denned in the event ((' — > (aa)(b(3), where ( and (' have the same masses 
m^, a and b are visible objects, and a and (3 are invisible particles with the same mass M. 
In such a event, m T2 variable is defined as follows, 



m T 2{PT, Pt> F^t; M test) = a min 



max|m T (p^,^;M tes t),m T (p^,^;M tes t)|j . (9) 
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FIG. 2: mp H ^ (vertical) vs. fnp H (horizontal) distributions for T_T_ (the left figure) and for tt 
(the right figure). The units of axes are GeV. 
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FIG. 3: Apt/pt,h and AR for pn 1 - 

Here, M tcst is an arbitrary chosen test mass and the transverse mass my is defined as 
follows, 

m 2 T (p a T , M tcst ) = m 2 a + M t 2 cst + 2 [E^ - p^] . (10) 

It is important that the following condition is satisfied in the case M test = M: 

m T2 (M) < m c . (11) 

Thus can be extracted with measuring the upper endpoint of the m T2 (M) distribution 
(m™| x (M)) in the case that the true M is known. 

In the case the true M is not known, we can calculate rxiTi (M test ) for an arbitrary test 
mass Mtest. For each test mass M tcst , we can measure m™2 X (M test ). The end point is 
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expressed in terms of the following equation for the case that the masses of visible systems 



are the same (m V is) and there is neither initial nor final state radiation 44|, 



mj + ml-M 2 fm 2 + m 2 --M 2 \ 
m?r( M») = + ^ (_S_^ ) + Ml,, - (12) 

The endpoint contains the information on a combination of the relevant masses and 
M. For a particle which undergoes more complicated decay process, m v i S can be various 
values. Therefore one can obtain more than two independent information on the masses. 
Practically, one can extract the true mass M from a kink structure of the rriT2 endpoint as 
a function of M test [44 1. 



For our case, visible particles are two top quarks, and invisible particles are two Aff's. 
It is, therefore, not possible to determine top partner mass itself from the kink method 
because m vis is always m t and not to be various values. We show the m T2 { M test ) — M test 
2-dimensional scattering plot in the Fig HI The dashed line is the line defined by eq. ( fl2l) 
substituted with the nominal values. The test mass dependence of the endpoint m™2 X (M test ) 
is well described by eq. (JT2l) and no detectable kink structure can be seen. Eventually, we 
can measure only a combination of masses; 



+ mf — m\ 



(13) 



2rriT_ 

In the case that a system of pair-produced particles has a net transverse momentum, the 
m™| x (Mtest) changes greater than the eq. f|T2|) for all M tes t but for M tes t = mi fl , therefore 



a kink structure might be seen [45] . However, T_T_ system generally has a small net 
\pt\ ~ 0(100) GeV in average, a kink structure is not seen in our case even at parton level. 
Then we cannot measure the m-r_ unless the itla h is known for our case. The itla h may 
be determined if productions of the other T-odd particles are observed. Alternatively, if we 
assume the thermal relic density of An is consistent with the dark matter density in our 
universe, is related to the Higgs mass so that it is determined with two fold ambiguities 
46]. 



Next, we show that toji can be measured using reconstructed tops at jet level assum- 
ing m>A H is known. Ideally we may regard hemisphere momenta as top momenta if both 
hemisphere masses satisfy the condition 150 GeV < rrip H < 190 GeV. There is, however, 
not enough number of events left under these cuts as mentioned in the previous subsection. 
Therefore in the following, we apply the cuts that one of the hemisphere masses mp H satisfies 
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FIG. 4: left) (Attest) — Attest 2dim. scattering plot at the parton level, the dashed line is the 
line defined by the theoretical m™| x (M tes t) description with no transverse momentum limit, right) 
The tbt2 distributions for the nominal mass rriA H = 151.8 GeV for the "or" (solid line). The 
dashed line shows the contribution from the events survived after the "both" cut. The endpoint is 
801.0 ± 9.4 GeV, and this value is consistent with = 800 GeV. 

150 GeV< mp H <190 GeV while the other satisfy 50 GeV < mp H < 190 GeV (the column 
"or" in the table) and regard the hemisphere momenta as top momenta. 

The endpoint of the rriT2 distribution does not change under the relaxed cut "or" . because 
is an increasing function of visible masses. Additional sources of missing momentum 
(such as neutrinos) do not affect the endpoint either. It is easy to understand this as follows, 
a system of a LTP and the other sources of missing momentum can be regarded as an invisible 
pseudo-particle. The invisible pseudo-particle's invariant mass is always larger than the LTP 
mass (m invisib i c = a/(^ltp + Pother) 2 > ™>Ah)> and m T2 ("^invisible) is nevertheless smaller than 
because the system comes from QC,' pah production. On the other hand, m^^ftest) is 
a monotonically increasing function of M tes t, therefore, mnlm^) < m^2 (^invisible) < "\ is 
satisfied. 

The rriT2 distributions for the nominal value vtia h — 151.8 GeV are shown in Fig. Upright). 
The distribution after the "or" cut for the signal events are shown in the solid line. We fit 
the distribution near the endpoint by a linear function and obtain m™2 X = 801.0 ± 9.4 GeV. 
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This value is consistent with the nominal value rriT_ = 800 GeV. 3 The fact supports validity 
of the relaxed cut ("or" cut) in determination of the rriT2 endpoint. The contribution from 
the events survived after "both" cut is shown in the dashed line and there are a few statistics. 

The distribution for the SM backgrounds is also shown in a dark histogram, and they 
have lower values. After imposing the cut m^ 2 > 350GeV, the SM background is 
significantly reduced but the signal is not reduced as in Table [H Moreover, after imposing 
the cut rriT2 > 500GeV, the SM background becomes negligible. Therefore we can neglect 
them to fit the endpoint ~ 800GeV in the present case. 

In a case that is lighter, for example = 500GeV, a top quark momentum arising 
from a top partner is approximately 200GeV [jri^ H = 150GeV is assumed). It is boosted 



enough and the decay products distribute within AR < 1.5 [471 ]. then the hemisphere ana' 
may work well. And the <r(T_T_) is ten times larger than the case of m-r_ = 800GeV 
Therefore it is possible to measure the endpoint by the method we proposed above 4 . 



ysi s 



25]. 



III. COMPARISON AMONG JET RECONSTRUCTION ALGORITHMS 
A. Jet reconstruction algorithms 

We now study a dependence of the signal distributions on jet reconstructing algorithms. 
The reason to study different jet reconstruction algorithms is as follows. Note that we need 
to study the jet system arising from a boosted top quark. In the rest frame of a top partner, 
a top quark momentum arising from top partner decay is expressed as 

p t =—- x l- 2 - + — - 365 GeV, (14) 

2 y m4_ 

therefore, typical px of a top quark is above 300 GeV. The jet angle separation is of the order 
of 9 ~ ttit/pt- If decay products of a top are aligned in the direction of the top momentum, 
the angle is even smaller. It is important to choose the algorithm that gives the best result 
in such a situation. Four algorithms (Snowmass cone, kt, Cambridge, SISCone) are used in 
the following analyses. 

3 In Figf4l we do not include the effects of jet energy smearing. In Appendix I A 11 we discuss the point. 

4 For discovery, we can also use lepton channels. In Ref. [4a], it is found that T_ can be discovered in the 
case of AMta = »^t_ — wia h ~ 250GeV. 
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1. Cone algorithms 

We take two cone-type algorithms. The first one is "Snowmass cone", which is a simple 
algorithm implemented in AcerDETl .0. It defines a list of jets as follows, 

1. Find the particle % which has the maximum E in all particles, and take it as 
a seed. If Ei is less than some threshold E t ^ then the process is finished. 

2. Sum the four-momenta of the particles in the circle of whose center and radius 
are (rj seed , ^seed) and R respectively. Define the four- momentum as p con e and 
redefine it as a new seed. 

3. Repeat Step 2 until p CO ne is converged. a 

4. Remove the constituents of the cone from the particle list and repeat from 
Step 1. 

a AcerDET jet finding algorithm skip this iteration. 

In this paper, we take cell momenta as massless particle momenta. As one can easily see 
from the algorithm, the highest p? jet in a region takes all activities within R(< 0.4) even 
if there are sub-dominant activities nearby R ~ 0.4. As we will see later, jet-parton energy 
matching is worse than the other algorithms. 

The second one is SISCone j^. This algorithm is a seedless cone search algorithm. It 
defines a list of jets as follows, 
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1. Find all "stable" cones seedlessly and calculate four-momenta of these cones. 
Here, a "stable" cone is defined with a set of particles satisfying the following 
condition, 

d(pi,Pcone)<R 

d(pi,pj) is the distance in the r\ — plane. 

2. Remove cones which have less energy than some threshold E t ^ from the cone 
list. 

3. If there are overlapping cones, determine to split or merge according to the 
overlap parameter /. Namely if the fraction of overlapping activities of the 
two jets by the smaller jets is larger than / two jets are merged, otherwise 
split the overlapping activities into the two jets. And update the cone list. If 
there is a cone which is not overlapping with other cones, remove it from the 
cone list and add it to the jet list. 

4. Repeat Step 3 until there is no cone in the cone list. 

The reconstruction algorithm is infrared safe, because the reconstruction does not relay on 
the highest pt cell in a cell list. The number of reconstructed jets depends sensitively on 
the overlap parameter /. If we set / smaller, the algorithm tends to merge jets. For our 
choice / = 0.75, which is the default value, jets from a top quark tend to be merged and 
efficiency of resolving the three jets in a hemisphere is rather low compared with the other 
algorithms. 



2. Clustering algorithms 



The other category of jet finding a 



in this category, the kt algorithm 
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gorithms is a clustering algorithm. A typical algorithm 
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371 ] is defined as follows: 
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1. Work out the kt distance dij for each pair of particles with momentum hi, kj 
and diB for each particle i. 

d l3 ee min(4, fc* d iB = kl R% ee (A V )l + (A0)J. (16) 

2. Find the minimum (i m i n of all the dij, diB- If ^min is a d^-, merge the particles 
% and j into a single particle by summing their four-momenta. If the d m { n is a 

then regard the particle i as a final jet and remove it from the list. 



3. Repeat from Step 1 until no particles are left. 



Cambridge algorithm 38|, |39( is similar to the kt algorithm but definition of d^ and di B is 
modified as follows: 









a %3 - R2 , 


d iB = 1- 


(17) 



B. Hemisphere invariant mass distributions 

In Table HI1 we show the numbers of signal events for the four jet algorithms after the 
same hemisphere mass cuts as in Table HI AcerDET has an option to rescale jet energy, the 
results for the Snowmass cone algorithm are given with jet calibration. The scale factor 
is determined so that an invariant mass distribution of the two jets from W has the peak 
consistent with mw 



For T_T_ events 



Snowmass cone (calibrated) 
kt 

Cambridge 
SISCone 



375 
420 
404 
425 



m t 



442 
415 
396 
385 



both 



143 
135 
130 
137 



363 
411 
398 
396 



m,T2 endpoint 



795 GeV 
797 GeV 
801 GeV 

796 GeV 



TABLE II: Summary of the cuts for various jet finding algorithms. The cuts are the same as in 
Table E 

We obtained 375, 420, 404, 425 events after the mp Hi ~ m t cut for the Snowmass cone, 
kt, Cambridge, SISCone respectively 5 . 



As we mentioned already, the Snowmass cone algorithm in AcerDET ignores jet invariant masses, therefore 
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FIG. 5: The distributions of mp H ^ for the T_T_ events for a) Snowmass cone 
(AcerDET) , b) kt, c) Cambridge and d) SISCone. 

The distributions of mp Hi for four algorithms are shown in Fig. [5j The shaded regions 
denote the region satisfying 150 GeV < mp H < 190 GeV. The kt, Cambridge and SISCone 
show nice resolutions in top mass. The peak for the Snowmass cone algorithm is dull 
and has broad tail. This is because AcerDET takes massless jets. This is rather an artificial 
difference as it is straight forward to define non- zero jet masses from calorimeter information. 
If such jet definitions are feasible at the LHC environment, the reconstruction efficiency may 
be increased significantly although we have not simulate the effect of mis-measurement of 
calorimeter energy. We regard the efficiency in the Snowmass cone as a conservative estimate. 
Fortunately, the endpoints of distributions are rather insensitive to the reconstruction 
algorithm. We find that they are 795 GeV (Snowmass cone), 794 GeV (kt), 801 GeV 
(Cambridge), 796 GeV (SISCone) with the statistical errors of the order of 10 GeV. The 



jet energy calibration should increase a jet energy to compensate the missing jet mass. The calibration also 
compensate average energy of particles that fall outside jet cones. For our case, several jets go collinear, 
and the particles outside the cone often fall into other jet cones, leading overestimate of the jet energies. 
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difference among the algorithms is not essential at this point. 
C. Parton-jet matching 




-0.2 -0.1 0.1 0.2 -0.2 -0.1 0.1 0.2 

Snowtnass kt 




-0.2 -0.1 0.1 0.2 -0.2 -0.1 0.1 0.2 



Cambridge SISCoiw 

FIG. 6: Ap<r/pT,Hi distributions of hemisphere H% for T_T_ events. The events with 150 GeV< 
mp Hi < 190 GeV are selected. 

Fig. [6] shows deviation of a hemisphere momentum from a true top parton momentum 
APt/Pt,Hi for the four algorithms. We selected the signal events with 150 GeV< mp H < 
190 GeV. The APt/Pt,^ is mainly distributed within ±5% region for all algorithms. The 
AR for top parton and hemisphere momentum (= a/ (tjh — Vtop) 2 + (4>h — 0to P ) 2 ) is mainly 
distributed less than 0.03 for all algorithms. These agreements justify regarding a hemisphere 
momentum as a top momentum. The peak position is larger than by approximately 2% for 
the Snowmass cone algorithm, because jet energy calibration by AcerDET leads over-estimate 
of the jet energy for collimated jets. The jet energy calibration compensates the activities 
outside the jet cone, but for the collimated jets they are taken into account by the other 
jets. For the other algorithms, the peak position is less than due to semi-leptonic decays 
of b- quark. 

We now look into the matching between a 6-parton and a 6-jet in a selected hemisphere. 
The selected hemispheres satisfy the following conditions: 1) there are only three jets in 
the hemisphere and 150 GeV< rrijjj < 190 GeV, and 2) at least one jet pair satisfies 
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Cambridge 



FIG. 7: ApT,b-jet/PT,b-paxt distributions for the selected hemispheres (150 GeV< mp H < 190 GeV) 
for the T_T_ events. 
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mw\ < 20 GeV and the other jet is 6-tagged. Here, we define a jet with AR b < 0.2 



as a 6-jet (AR b = a/ (rj iet - rj b - paxt ) 2 + (0 jct - 0b_ part ) 2 ). In Fig. U\ the Ap T)b _ ]ct /p T ^ pa , rt 
distributions are shown. The large tails found for Ap Tb _j et /p Tb _ paiTt < come from a 
6-parton decaying semi-leptonically. In addition, the distribution for the Snowmass cone 
algorithm shows a tail for Ap Tb _j et /p Ttb _ paiTt > 0. This tail arises because a locally highest 
Pt jet takes over all energy in R = 0.4 cone because the jet finding algorithm starts from 
the highest pt clusters. This feature cannot be improved with minor modifications of the 
algorithm. For the SISCone algorithm, the number of the three jet events is significantly 
small compared with the others because the algorithm actively merges overlapping jets. We 
will discuss this point in the next subsection. The algorithm is therefore not suitable for our 
analysis in the next section, in which we study the top spin dependence of 6-jet distributions. 



D. Number of Jets distribution 



In this subsection we compare clustering algorithms with the SISCone in terms of the 
number of jets in a hemisphere. The numbers of jets inside a hemisphere with 150 GeV< 
mp H < 190 GeV are shown in the Table HTT1 and HVl Since kt algorithm behaves similar to 
Cambridge algorithm, only those for the Cambridge and SISCone are shown. The parameter 
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R for the clustering algorithms and for the SISCone have different meanings and SISCone 
has additional parameter / as explained in Sec lHI Al We investigate the distribution of 
the number of jets varying these parameters (R = 0.2, 0.3, 0.4 and / = 0.5, 0.75, 0.9 for 
SISCone). 

In order to study the top decay distribution, it is better to choose the parameters which 
give higher 3-jets acceptance. With the Cambridge algorithm (R = 0.4) 511 events are 
classified into a group of 3-jets events, while with the SISCone (R = 0.4, / = 0.75) 324 events 
are classified into it although the total numbers of hemispheres with 150 GeV < mp H < 190 
GeV are approximately the same. We can see (R, f) = (0.3, 0.9) or (R, f) = (0.2, 0.5) 
are optimal to enhance the number of 3-jets events with the SISCone for our model point. 
The distribution of numbers of jet at those parameters are similar to that with Cambridge 
(R = 0.4). For such a small R, however, some activities are missed outside a jet cone 
leading worse parton-jet matching 6 . Moreover, we found Cambridge with R = 0.3 also gives 
higher acceptance for 3-jets events than with R = 0.4. We do not find the parameter which 
improves the results for SISCone over Cambridge by changing R and /, therefore, we use 
clustering algorithms for the further analysis. 

Appropriate R should be used depending on top to protect unnecessar y m erging. Sub- 
jet analysis based on clustering algorithms might be useful in such a case [49j. We do not 
discuss these points any more because it is beyond the scope of this paper. 



R = 0.4 


Cambridge 


SISCone 


/ = 0.5 


/ = 0.75 


/ = 0.9 


1 jet 


9 


103 


58 


41 


2 jet 


244 


429 


413 


362 


3 jet 


511 


295 


324 


362 


4 jet or more 


36 


5 


15 


23 


total 


800 


832 


810 


788 



TABLE III: Summary of the number of jets in a hemisphere with mp H ~ m% for R = 0.4. 



6 For such a small R, detector granularity might not be enough to resolve the jet 
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Cambridge 


SISCone R = 


0.3 




Cambridge 


SISCone R = 


0.2 




J? = 0.3 


f = 0-5 


/ = 0.75 


/ = 


0.9 


R = 0.2 


f = 0.5 


/ = 0.75 


/ = 0.9 


1 jet 


2 


51 


16 




12 


1 


10 


2 


2 


2 jet 


138 


335 


294 




235 


47 


135 


106 


84 


3 jet 


574 


419 


430 




462 


467 


477 


440 


430 


4 jet or more 


85 


36 


55 




71 


136 


107 


122 


115 


total 


799 


841 


795 


780 


651 


729 


670 


631 



TABLE IV: Summary of the number of jets in a hemisphere with mp H ~ mt for R = 0.3 and 0.2. 
E. Effects of Underlying Events to the reconstruction 




mHl mHl mHl mHl 



Without underlying events With underlying events 

FIG. 8: The mp^ distributions in Cambridge (solid) and kt (dashed) algorithms without (left) 
and with (right) underlying events. We take R = 0.4 ~ 0.7. 

So far we have discussed the event distributions without underlying events. Underlying 
events come from the soft parton interactions which occur with a hard collision, and whose 
nature at the LHC has large theoretical uncertainty. The top reconstruction efficiency 
may become worse with them, because the number of hit cells significantly increases with 
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underlying events. We have generated the signal events with underlying events and multiple 



scattering using HERWIG6 . 5 + JIMMY 4l|. In Fig. [8] (left), we show the distributions of mp H ^ 
for the kt (dashed) and Cambridge (solid) algorithms without underlying events for R = 0.4, 
0.5, 0.6, 0.7. The event selection cuts are the same as section III Bl We can see that the 
locations of the peaks increase as increasing R for both the kt and Cambridge algorithms. 
The shapes of distributions are similar for all R as top quarks are boosted enough so that 
the decay products are isolated from the other activities. However the kt algorithm tends 
to give higher invariant mass than the Cambridge algorithm. 

In Fig. [8] (right), we show the same distributions with underlying events. We can see 
that the position of the peak for the kt algorithm is significantly larger than that for the 
Cambridge algorithm in all R values. Even for R = 0.4, the peak position is larger than 
175 GeV for the kt algorithm. The situation becomes worse for the kt algorithm as R 
increases. The reconstruction efficiency is reduced significantly for R > 0.6. This is because 
the kt algorithm over-collects soft activities which are far from the jet direction (large Rij) 
due to the factor min(fc^,fc^) in the definition of the distance in Eq. (116j) . which is known 
as splash-in effects. On the other hand, the Cambridge distance measure does not have the 
factor, therefore, it is not too sensitive to the existence of the underlying events. The effect 
of the underlying events can be safely neglected for R ~ 0.4. Hence we take the Cambridge 
algorithm in Sec. [Til and IIV1 

IV. TOP POLARIZATION EFFECTS 

In this section, we consider top polarization effects. In the Littlest Higgs model with T- 



parity, the Lagrangian relevant to a top partner decay is written as follows [15|, |20j, |28|, |46|, [50] . 



2g' _ 

C = i— cos 9 (sin (3Pl + sin a Piijt, (18) 





where, 



m t v mfv 

sina~ -, ship ~ — 5 — -. (19) 

m T _f m l T J 



A top partner T_ decays dominantly into a top with h t = +1/2 if 3> vtia h + n^t and 
sin/3 <C sin a, where h t is defined as the top helicity. It is the case in our model point, since 
sin/3 ~ rn t /rriT_ sina ~ 0.22 sina. The amplitude is calculated in the Appendix IA 21 and 
we find V = [N(h t = 1/2) - N(h t = -l/2)]/[N(h t = 1/2) + N(h t = -1/2)] ~ 0.85. 

22 



To simulate the top polarization effect we need to follow the decay cascade till the partons 
arising from top decay. Instead of generating pp — > T_T_ — > bbW(-^> qq)W(-^ qq)AnAH 
using COMPHEP, we generate stop pair production pp — > tt* followed by t — > t\\ at a MSSM 
model point using HERWIG 7 . We take the MSSM parameter that the other sparticles are 
heavy and the decay vertex of t\ is approximately proportional to t*X\tR, so that HERWIG 
generates stop pair efficiently and they decay into approximately completely polarized top 
quark (V = 0.996 for our model point). 

HERWIG has an option to switch off polarization effects. For the MSSM point, we do not 
find any distinguishable difference between the 777^2 distributions with/without polarization 
effects. Therefore the results shown in SecJTT] may be valid even for the LHT because the 
spin correlation effects are not large. 



Decay distributions of the top quark contains information on the interaction vertex 
The amplitudes for h t = ±1/2 are expressed as follows, 

jgrcosfe* (h t ,X w ,h b ) = (+,0,-), 
-v^sinfe 2 ^ (+,-,-), 



5l| 



M ~ \/2m t E b x < 



(20) 



-^sinf (-,0, 



>/2cos§e* (-,-,-)• 

The amplitudes for the other helicity combinations vanish. Here, E b = (mf — m^ v )/2m t is 
the energy of the bottom quark in the rest frame of the decaying top. And 9 and are the 
polar and azimuthal angles of the momentum of W boson defined in the rest frame of the 
decaying top. The 9 is measured from the top momentum direction and the Xw is a helicity 
of W. 

In Figj9] (left), we show the parton level top decay distribution as the function of cos 9 tb for 
the polarized and non-polarized cases, where 9 tb is the angle between the 6-quark momentum 
and the t-quark boost direction at the rest flame of the t-quark. We find evident difference 
between the two cases. In the region around cos 9 tb ~ 1, the emission of b quark is suppressed 
because the amplitude proportional to cos(6*/2) is suppressed for h t = +1/2. 



7 Note that the spin of the intermediate particle are different between these two processes. Especially, 
we expect spin correlation between T_T_ and their decay products, which does not exist for t\t\. The 
correlation in principle appear in momentum distribution of t and t. However, the effect is rather small 
because T_ is non-relativistic and also the system does not have enough kinematical constraints. 
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polalized ■ 
non-polalized ■ 



FIG. 9: left) The polar angle (cos 0^6 = cos(7r — 9)) distributions of bottom quarks from polarized 
(solid) and non-polarized (dashed) top decays, right) The polar angle (cos#^*) distributions of the 
6-tagged jets. 

To see this at jet level, we study distributions of jets that are consistent to top decay 
products. Some parton configuration is difficult to resolve at jet level. We analyse only 
the hemispheres with 150 GeV< mp H < 190 GeV. We require that two of the three jets 
are consistent with those coming from W, that is \rrijj — mw\ < 20 GeV, and the other 
jet is 6-tagged. Here, we regard a jet as a 6-jet if the direction of the momentum is in a 
AR = 0.2 cone centered at a bottom parton momentum with > 20 GeV. The analysis of 
the hemisphere which consists of 2 jets is given in Appendix IA 31 

In FiglH] (right) we show the distribution of the angle between the b-jet momentum 
and the reconstructed top momentum Off. For the plot, we selected only the events with 
m<T2 > 500 GeV. Under the cut, Standard Model backgrounds are negligible as seen in SecJITI 
We use the measured hemisphere momentum to go back to the rest frame of the jet system. 
There is a distinguishable difference between polarized and non-polarized distributions. The 
ratio n(cos 9lf < O)/n(cos0jf > 0) = 2.08 for the polarized case, while it is 1.16 for the 
unpolarized one. 

A polarized top quark decays into a polarized W (\w = 0,—). Decay distribution of 
polarized W (W — > 2j) can be calculated and the amplitudes are written as follows, 



M- oc 



M oc 



cos 9* 



-icj>* 



sin i 



V2 



1 + cos#* ir 
M+ oc e l<p 



(21) 
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Here, 6* and 0* are the polar and the azimuth angles of the momentum of one of the jets 
from a W decay to the W momentum direction. These momenta are defined at the rest 
frame of the W. A longitudinally polarized W (X\y = 0) tends to decay transversely. On 
the other hand a transversely polarized W (Xw = ±) tends to decay along a direction of the 
W momentum. 

These differences may appear in the jet pt asymmetry A, which defined as follows, 

aJpILZIiA. (22) 

PTl + PT2 

Jets from a W(Aw = 0) decay tend to have A ~ while those from a W(X = ±1) give 
larger A. This can be seen in FigJTU] (left). In this plot the A pSuT t distributions are shown for 
the events with cosfljf < 0. The ratio N(W{X W = 0))/N{W{X w = ±1)) with cos 0*6 < is 
larger for polarized tops. Therefore, A for polarized top is distributed more around than 
for non-polarized top. Fig. [TD] (right) shows A distributions at jet level for cos9 3 t f < 0. 
Unfortunately, It is difficult to see the differences only by the shape of the distribution due 
to the limited statistics. 




FIG. 10: The distributions of the pt asymmetry A of the jet pairs from W decays for polarized and 
non-polarized cases. Parton level (left) and jet level (right) distributions are shown. The difference 
of the distributions can be seen even at jet level. 



V. CONCLUSION 

In this paper we have studied reconstruction of top quarks arising from T_T_ produc- 
tions and its subsequent decay into a top and a stable gauge partner A H in the LHT. We 
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demonstrate the reconstruction of the top quarks through finding collinear jets whose invari- 
ant mass is consistent with m t using hemisphere analysis. Main SM background processes 
are tt+jets, Z+jets and W^+jets productions, which can be reduced by imposing the cut on 
hemisphere momenta and vtiti variable. 

We also investigate the dependence on jet reconstructing algorithms. The cone algorithm 
used in the previous study [28| is not optimal for the process. A top from T_ is boosted, 
while the algorithm is designed so that the highest pt jets take all activities near the jet, 
mis-estimating the energy and the direction of the jets. An infrared safe version of the cone 
algorithms (SISCone) also has some disadvantage for our case, because they tend to merge 
overlapping jets. We also study distributions with modern clustering algorithms (kt and 
Cambridge), which in general give better results than the cone algorithms. 

We also study effect of underlying events, and find the known tendency that the kt 
algorithm overestimates jet energies caused by collecting far and soft activities. Whereas 
the reconstruction efficiency in the Cambridge algorithm is not affected if R ~ 0.4. 

We also discuss top polarization effects. A top quark arising from a top partner decay 
is naturally polarized. This can be studied through looking into a distribution of the 6-jet 
from t decay especially the angle to a reconstructed top momentum in the rest frame of 
the top. We find that difference of the distributions between polarized and non-polarized 
top is still at detectable level with 6-tagging for reasonable integrated luminosity (J C =50 
fb _1 for rrtT_ = 800GeV and rriA H = 150GeV). These analyses are demonstrated using the 
Cambridge algorithm, which shows good 6-jet and 6-parton matching. 

In many new physics scenarios, boosted gauge bosons and top quarks are produced at 
a significant rate. Our study shows that choosing a right jet reconstruction algorithm or 
studying the dependence on them is important to reveal the physics behind the signal. 
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APPENDIX A: APPENDIX 
1. Jet smearing 

In this paper, we show the distributions without jet energy smearing. In AcerDET, there is 
an option to smear jet and missing transverse energies. The smearing is introduced for each 
jet energy in the Snowmass cone algorithm and for a sum of the total transverse momentum 
(the missing transverse momentum) rather than for each calorimeter cell. We do not try to 
include smearing effects for the other three jet reconstruction algorithms (kt, Cambridge, 
SISCone) in this paper, because cells that jets consist of depend on the jet reconstruction 
algorithms, therefore comparison of smearing effects under the same condition is not easy. 

To obtain a rough idea on the signal and background distributions with smearing, we 
show m T2 distributions with/without smearing and jet energy calibration in the Snowmass 
cone algorithm. The distribution near the end point is not significantly changed. The 
effect is rather small because we take the smearing based on ATLAS detector performance 
50x v / ^%. which is less than 10 % for jet with p T > 30 GeV. 



FIG. 11: Distributions of m.T2 for the Snowmass cone algorithm without/with smearing (left /right 
figure respectively). The dashed line shows for the events with 150 GeV < mp Hi < 190 GeV for 
both H{. The solid line shows for the events with 50 GeV < mp Hi < 190 GeV. The endpoints are 
785.9±6.6/810.2±25 GeV for nominal m All = 151.8 GeV . (m T _ = 800.2 GeV). Calibrated jets are 
used for both figures. Dark histogram is the background distribution in the case that jet energy 
smearing is on. We do not show the background distribution without jet energy smearing in the 
right figure. 
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Smearing on 



Smearing off 
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2. Top polarization 



In this paper we took mass parameters m T - = 800.2 GeV, vtla h = 151.8 GeV, and 
m t = 175. GeV. The Lagrangian relevant to our study is 

2q' - 

C = i-?- cos9 H T-4 H (sm (3P L + smaP R )t. (Al) 
5 



Here, a and (5 is approximately expressed in terms of 



,2, 



sina~ -, sin (3 ~ 2 , (A2) 

m T _f m 2 T J 

therefore sin/3 ~ 0.22 sin a at our model point. 

The amplitude of a top partner decay into a top with helicities h t , hx, \a can be calculated 

as follows, 

M ht , hT ,x A i< tA H \i4 H {s p P L + s a P R )T_|T_ > 

= - ie h A ^VA H ]rn AH )u ht (-p t ]rn t )^(spP L + s a P R )u hT (p T _; m T ) 



= -i e i(h T -h t +x A )4> ^/^T 



-^H^^WC^^i + 2fc*^4i)] , (A3) 
where at the rest frame of T_, 

Pt — (ptsin0cos0,p t sm#sin0,p t cos#) = — Pa h , 

E, = ^ + A, £. = ^-A, (A4) 

A = =Ll^, (A5) 
zmr 



m T _ m 2 t +m A (m 2 t - m\ ) 2 

Pt = PA H = ^-\\l-2 -2 + -4 , (A6) 



and, 



A* ± = v 7 ^ ± ™*, (A7) 
1a = E A[1 /m AH , (3 a =p Ah /E Ah . (A8) 

i 

The h2 (6>) is the Wigner's function, 



(A9) 
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The square of the amplitudes are therefore expressed for each helicity eigenstate as follows; 
2 sin 2 f(^AL + ^M' + ) 2 , 

a+££ j^t _|_ Sq-Sff y^f J _|_ ^ go+gff ^t _|_ s a -sp ^ 

+ ^M' + )] 2 , (+, 



2cos 2 f( 



7 icos 2 |[^(^ 



( + ,+ 
12 ( + : 



7 2 sin 2 f 

|7W| 2 ~ 2m T _ hsin 2 §(^±^4*_ - ^ 
2cos 2 f(^M*_ - 



7isin 2 f[/5 A (: 



7 2 cos 2 |[^(^ 



I 2 , 



^ + )] : 



(■ 



-,+), 
,-), 

+,o), 



- 0), 
-,+)(AlO) 



-,+,-), 

-,+,o), 
-,o), 



(the others). 

Polarization of top quarks from top partner decays is given as follows, 



v = r(T-^t (+) A g )-r(T-^t(_)A H ) 

= (si- g 2 )M(2 + 7i + 7l/gj) 

~ (* 2 + 4)( 2 + 7i + W - 2^(2 + 7l - 

We obtain P ~ 0.85 for our model parameter. 



(AH) 



3. Jets Configurations in Hemispheres 

For hemispheres consistent with a top mass, some parton configurations are difficult to 
resolve at jet level. We categorize hemispheres with 150 GeV< mp H < 190 GeV into the 
following four groups: 

I. Only one jet in a hemisphere. 

II. Only two jets in a hemisphere. One of the jets has mass consistent with W(m w — 20 GeV 

< rrij < m w + 20 GeV) and the other is 6-tagged. (W decay products are merged into 
a jet and 6-jet is isolated) 
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FIG. 12: The distributions of the 6-partons in hemispheres that the total hemisphere mass is consis- 
tent with wit as the function of 0^. The right(left) figure is for (non-)polarized tops. Distributions 
are shown for the four different groups. I) one jets in the hemisphere (dashed thin line). II) two 
jets in the hemisphere and one jet has mass consistent with W (dashed thick line). Ill) two jets 
in the hemisphere but none of jets has mass consistent with W (solid thin line). IV) three jets in 
the hemisphere and mass of a pair of two jets is consistent with W (solid thick line). 

III. Only two jets in a hemisphere. None of the jets has mass consistent with W mass and 
at least one of jets is 6-tagged. (One of the partons from W decay and a 6-parton are 
merged into a jet, and the other from W is isolated. ) 

IV. Only three jets in a hemisphere. Invariant mass of two jets is consistent with W and 
the other jet is 6-tagged. 

Here, we regard a jet as a 6-jet if the direction of the momentum is in a AR = 0.2 cone 
centered at a bottom parton momentum with p? > 20 GeV. 

In Fig. [121 we plot 6tb distribution of jets for each group, where 6 t b is the angle between 
a 6-parton and a t-parton momenta at the rest frame of the t-parton. The left figure is for 
non-polarized top quarks and the right figure is for polarized top quarks. 

For type IV events, the number of the events where a 6-parton goes in the forward 
direction is strongly suppressed for polarized events compared with non-polarized events, 
while the events for cos 9 t b < is significantly enhanced. It is consistent with parton level 
distributions. Note that the events near cos 9 t b ~ — 1 cannot be accepted as type IV for 
both polarized and non-polarized cases because pt of the 6-parton is too small, while the 
region is the most sensitive to polarization effects. The other distributions do not show clear 
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